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Abstract 

We consider the evolution of N bosons interacting with a repulsive short range pair potential 
in three dimensions. The potential is scaled according to the Gross-Pitaevskii scaling, i.e. it is 
given by N 2 V(N(xi — Xj)). We monitor the behavior of the solution to the iV-particle Schrodinger 
equation in a spatial window where two particles are close to each other. We prove that within this 
window a short scale interparticle structure emerges dynamically. The local correlation between 
the particles is given by the two-body zero energy scattering mode. This is the characteristic 
structure that was expected to form within a very short initial time layer and to persist for all 
later times, on the basis of the validity of the Gross-Pitaevskii equation for the evolution of the 
Bose-Einstein condensate. The zero energy scattering mode emerges after an initial time layer 
where all higher energy modes disperse out of the spatial window. We can prove the persistence 
of this structure up to sufficiently small times before three-particle correlations could develop. 

Keywords: Bose-Einstein condensate, Gross-Pitaevskii equation, dispersive bounds. 
AMS classification numbers: 81U05, 81U30, 81V70, 82C10. 

1 Introduction and main result 

We consider a three-dimensional system of N indistinguishable spinless bosons coupled with a pair- 
wise repulsive interaction Vn- The Hamiltonian of this system is 

N 

h n = y,(-^)+ Yl Wxi-xj) (i.i) 

acting on L^ ym (M. 3N , dx) , the symmetric subspace of the tensor product of N copies of the one- 
particle space L 2 (IR 3 ). Here x = (xi,X2, ■ ■ ■ , xjy) G M. 3N denotes the position of the N particles. We 
assume that the interaction potential Vn scales with N according to the scaling introduced by Lieb, 
Seiringer and Yngvason in [BJ; that is, we fix a repulsive potential V and we rescale it by defining 

V N {x) := N 2 V{Nx). (1.2) 
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We assume that the unsealed potential V : R — ► M is non-negative, smooth, spherically symmetric, 

v sym( 



and compactly supported. The wave function of the system at time t is denoted by ^Nt S 



with H^at^I^ = 1. It evolves according to the Schrodinger equation 

idt* N ,t = H N V N , t (1.3) 

with a given initial condition ^^,4=0 = ^N- We will be interested in the evolution of initial data 
exhibiting complete Bose-Einstein condensation. 

For a given wave function $jy £ L 2 ym (M 37V ) we define the one-particle marginal 7^ = 7^ 
associated with tp n to be the positive trace-class operator on L 2 (M 3fc ) with kernel given by 

7^(x;x) := / ^ N (x, z 2 , ■ ■ ■ , z N )^ N (x' , z 2 , ■ ■ ■ , z N ) dz 2 ■ ■ ■ dz N . (1.4) 

Je3(JV-1) 

We say that a sequence {^tvItven of iV-body wave functions exhibits complete Bose-Einstein con- 
densation in the one-particle state (p G L 2 (M 3 ), ||v?||2 = 1 ; if 

^ ] ^\<p)(<p\ (1-5) 

in the trace-norm topology, as N — > 00. Here denotes the orthogonal projection operator onto 

(p (Dirac notation). 

On the level of the one-particle density matrix, the condition (jl.5)l is a signature that almost all 
particles (up to a fraction vanishing in the limit N — > 00) occupy the same one-particle state, de- 
scribed by the orbital <p (condensate wave function). Suppose now that a trapping external potential 
is added to the Hamiltonian, i.e., instead of (jl.ip consider the Hamiltonian 

TV 

H T P = ^(-Ai + [/(xO) + Yl V N (xi- Xj ) (1.6) 

i=l l^Kj^N 

with U(x) — > 00 as |x| — > 00. In [6j, Lieb, Seiringer and Yngvason proved that the ground state 
energy of divided by the number of particle N (the ground state energy per particle) converges, 
as N — > 00, to the minimum of the Gross-Pitaevskii energy functional 

8(<p) = I (|V^| 2 + C/M 2 + 47raM 4 ) (1.7) 

over all cp G L 2 (R 3 ) with \\1pW2 = 1- The coupling constant a is the scattering length of the unsealed 
potential V. Later, in [5j, Lieb and Seiringer also showed that the ground state of H^ p exhibits 
complete Bose-Einstein condensation into the minimizer of (11.71) . 

We recall that the scattering length of a potential V is defined as 

a : = lim |x|a;(x) = — / dxV(x) (l — uj(x)) , (1-8) 

M->00 8-7T J R 3 

where 1 — u is the unique non-negative solution of the zero-energy scattering equation 

(- A+ \V)(l-u) = (1.9) 

with boundary conditions u(x) — > as |x| — » 00. Several properties of u are collected in Appendix ICl 
It follows from the definition that the scattering length of the scaled potential Vjsr is ajy = a/N and 
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that the zero-energy scattering solution associated with — A + \Vn is the function 1 — u>n with 
wjv(i) = lo(Nx). In particular, ujn has a built-in structure at the scale A -1 . 

The time evolution of a condensate after removing the trap, i.e. setting U = 0, can be described 
by the solution of the Gross-Pitaevskii equation 

idtft = -A<Pt + 8na\(p t \ 2 (pt (1-1°) 

with initial condition tyt\ t= Q = </?• 

The first rigorous derivation of the Gross-Pitaevskii equation (jl.lOp from many-body Schrodinger 
dynamics (|1.3|) has been obtained in [2] under the condition that the unsealed interaction potential 
V is sufficiently small. This result was later extended to large interaction potentials in [3]. More 
precisely, it has been shown that if the family of wave functions {^atIatsn has finite energy per particle 
(in the sense that {^NjH^^n) ^ CN), and if it exhibits complete Bose-Einstein condensation in 
the one-particle state ip in the sense (|1.5p . then, at any time t ^ 0, VPjv,t = e~ tHNt ^iq still exhibits 
complete Bose-Einstein condensation. Moreover, the condensate wave function ipt at later times is 
determined by the solution of the Gross-Pitaevskii equation (jl.lOp with initial data ft=o = </>• 



At first sight, the condition (|1.5j) may indicate that wave functions describing condensates are 
very close to be factorized (for = <p® N , (|1.5p holds as an equality). However, the structure of the 
evolved wave function *$>N,t is much more complicated. In fact, it turns out that VPjv,t is characterized 
by a short scale correlation structure which plays an important role in the derivation of the Gross- 
Pitaevskii equation. If this short-scale structure were lacking along the evolution and the wave 
function Vl/jv,t were essentially constant in the relative coordinates X{ — Xj whenever |xj — Xj\ < i, 
then the condensate wave function ipt would solve the equation (jl.lOj) with coupling constant b (this 
is the first Born approximation of 87ra); we discuss this problem in more details in Appendix [El 

If the initial state ^jv has a short-scale structure characterized by 1— uj]\f(xi — xj) for \xi— Xj\ < jj, 
for example, it is the ground state of the system before removing the traps, then this structure seems 
to persist along the time evolution. The condensate wave function (pt lives on order one scales and it 
changes with time, while the correlation structure on the scale A -1 is conserved. However, Theorem 
2.2 of [2J states that even if the initial state does not have a built-in short-scale structure, the time 
evolution of the orbital wave function is still given by (jl.lOp . This indicates that the characteristic 
short-scale structure not only persists but dynamically emerges within a very short initial time layer. 

The main result of this paper is a description of the dynamical formation of this short-scale 
structure. Consistently with the underlying Gross-Pitaevskii equation, which is derived at any fixed 
t in the limit A —* oo, correlations in a large but finite system with A particles are expected to form 
in a short transient time (i.e., o(l) in A) and then to be preserved at any macroscopic time (i.e., 
times of order 1). 

In the previous works the presence of the local structure was captured by the -ff i^-energy estimate 



V-V 

1 - LU N {Xi - Xj) 
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2 dx ^ c(tftf,^* v ) (l.U) 



valid for every ^at G L 2 (M. 3N ) and every fixed indices i ^ j in {1, . . . A}. This inequality has been 
proven in [2] for small interaction potential (for large interaction potential, another a-priori bound 
on the solution ^N,t of the A-particle Schrodinger equation has been obtained in [3]). 

Under very general conditions on the trapping potential U, the r.h.s of (jl.lip is of order one for 
the ground state of H l ™ p . If ^jy did not carry a short-scale structure (1 — ujn{xi — Xj)) when 
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\xi — Xj\ ~ N 1 , say, it were essentially constant on this scale, then the main term on the l.h.s of 
(jl.lip would be of order 

I 7 r Vjv(x)| 2 dx~ / \V 2 u N (x)\ 2 dx ~ 0(N) . 

J\xi-xj\^l/N 1 _ UN(Xi - Xj) J\x\^l/N 

More generally, this shows that the short-distance behaviour of any for which the inequality 

{V N ,Hff* N )^CN 2 (1.12) 

holds, is asymptotically given by (1 — ujn(xi — %j)) in the regime where \xi — Xj\ ~ AT" -1 . Since 
is conserved along the dynamics, 

the same conclusion holds for ^jvt at any later time t > 0, thus for such initial data the local 
structure is preserved in time. 

The same argument cannot be used to detect the presence and the formation of the short-scale 
structure for states that do not satisfy the bound (|1.12p . This is the case for the completely uncor- 
rected initial state (p® N , with some (p £ H 2 (R 3 ), because (<p® N , H^Lp® N ) ~ N 3 (see Lemma EHJ). 
However, the validity of the Gross-Pitaevskii equation for the product initial state (Theorem 2.2 
of [2]) indicates that the short-scale structure, although initially not present, still forms within a 
very short time interval. The length of this transient time interval must vanish as iV —* oo since 
after the N — > oo limit, the Gross-Pitaevskii equation (jl.lOp is valid for any positive time t > 0. 
The proof of Theorem 2.2. of [2J still relied on (|1.12p after an energy cutoff iJjy ^ K N in the many- 
body Hilbert space. The energy cutoff artificially introduced the short-scale structure in the low 
energy regime. Although it was showed that, in the case k>1, the high energy regime ^ kN 
essentially does not influence the evolution of the condensate, the proof did not reveal whether the 
short-scale structure is indeed formed in the full wave function without the cutoff. The goal of the 
present analysis is to prove the dynamical formation of the expected pattern of correlations in the 
time evolution of an initially uncorrelated many-body state = l P® N ■ 

The above discussion suggests to compare &N,t with 

N 

] | (l - u N (xi - xj)) JJ<^(av) (1.13) 

l^i<j^N r=l 

after the transient time. Although (11.130 has the expected built-in short-scale structure in each 
relative variable, the true Af-body wave function at later time t > is much more complicated and 
it is beyond our reach to describe it precisely. The main reason is that it quickly develops higher 
order correlations as well; the typical distance between neighboring particles is of order A^ -1 / 3 , so 
within time t of order one, each particle collides with many of its neighbors. Moreover, the energy 
of the product initial state, 



is much bigger than the energy predicted by the Gross-Pitaevskii functional, by recalling that b > 8na. 
Thus there is an excess energy of order N in the system that remains unaccounted for. This excess 
energy must live on intermediate length scales for < 7 < 1 not to be detected either on the 
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local structure of order N 1 or on the order one scale of the condensate. Moreover, these excess 
modes must be sufficiently incoherent not to influence the evolution of the condensate. 

We cannot describe the evolution of these intermediate modes yet, we can only focus on the 
formation of the local structure (1 — ujsr(xi — Xj)). This is selected in a scattering process by having 



the smallest local energy. Our study is restricted to a spatial window where 



Xi 



< £ for some 



intermediate £ between the length scale N 1 of the expected local structure and the macroscopic 
order one scale of the system. If the initial datum is essentially constant in the variable far 
< £, then a short-scale structure emerges corresponding to the zero-energy mode 1 — ljn 



clS | X >i X j 

together with all higher energy modes; these quickly leave the ^-window while only the short-scale 
structure remains. We can thus monitor the formation of the local structure (1 — oon{xi — Xj)) 
in this window and prove its persistence for a while after its emergence. The time scale must be 
sufficiently large compared to the window size so that all higher energy modes disperse, but it has 
to be sufficiently small so that no three-particle correlations could develop yet. 

We provide a rigorous version of this picture in Theorem 11.11 We focus on the local structure in 
the relative variable X\ — x 2 ; by symmetry the same result holds for any pairs of relative variables 



Xi — xj . We define a cutoff function 



with 



if < r < 1 
if r > 2 . 



We consider the following time-dependent quantity 



6i(xi - x 2 ) 



*jv,t(x) 



1 - u N (xi - x 2 ) 



*jv(x) 



dx 



(1.14) 



(1.15) 



(1.16) 



where 1 — uj n is the zero energy scattering solution of — A + \ Vn-, and ^at(x) 
initial iV-body wave function. The cutoff scale £ is always assumed to satisfy 



II/li <p(?j) is the 



1 



(1.17) 



Here and throughout in the sequel we will make the convention that by An <C Bn one means that 
< (log N) k Ajy[ < CBn with a sufficiently large k and C. 

The decrease of J~t<[(t) has the natural interpretation of formation of a local structure at the scale 
iV -1 within the £- window. By controlling this quantity we demonstrate that the conjectured short- 
scale structure indeed forms within a short time of order N~ 2 and then it is preserved for longer 
times in a time window that is essentially iV -2 <C t <C iV"~( 2 ~io). Our main result is the following: 

Theorem 1.1. Let V : M 3 — » R be a non-negative, smooth, spherically symmetric, and compactly 
supported potential. Let Vn{x) := N 2 V{Nx) for N £ N. Let ip G L 2 (M 3 ) with \\ip\\ 2 = 1 and 

4 

m=0 

for some a > 3. Consider the Hamiltonian 

N 



1=1 
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acting on L 2 (M? N ), the initial datum ^n = <P® N , and its time evolution ^N,t '■= e ^w*^®-^ 
Consider the function J- Nit) defined in \1.16\) . Then 

r N{t) , crm ( 055*2* <^ + <£gl (1 „ gA ^ (1 . M ) 



N k Ni NH 
for all times t such that < t <C iV -1 . ^4s a consequence, 

F N (t) < ^tv(O) /or (A^) 4 < iV 2 t < Nio(N£)i . (1.21) 



i 



Remark 1. Eq. (|1.20j) does not look dimensionally correct (in our coordinates, time has the 
dimension of a length squared). The reason is that, to simplify the notation, we consider the length 
scale A characterizing the initial wave function tp, the radius R of the support of the (unsealed) 
potential V, and the scattering length a of V as dimensionless constants of order one. More generally, 
if, for A > we set 

(pW(x) = X~ 3/2 ifix/X) 
the bound Eq. (|1.20j) assumes the (dimensionally correct) form 

^ ()<c ^ (0) (^(^ + ^ Pog( „ W ) 

where the dimensionless constant C depends on the ratios R/X and a/A (a more refined analysis 
would also allow to compute the precise dependence on R and a). 

Remark 2. Since the interaction potential has a length scale 1/N, it will be convenient to introduce 
the length L = Nl, expressing the size of the window relative to the interaction range. The two 
particle scattering process takes place on a time scale 1/N 2 (see Section [2]), thus it is also natural 
to introduce the rescaled time T = N 2 t. The appearance of the combinations of N£ and N 2 t in the 
theorem is motivated by the fact that we actually describe a long time and large distance scattering 
process in terms of the rescaled variables T 3> 1, L ^ 1, in the regime where L 4 < T < N 1 / 10 ^ 2 
(in the spirit of Remark 1 above, we consider the regime (L/A) 4 <C (T/A 2 ) <C A^ 1 / 10 (L/A) 1 / 2 ). 

Remark 3. To be concrete, choosing £ = ■h, the short-scale structure given by 1— ujn is established 
for times 

w « f « ■ (1 ' 22 > 

The first inequality corresponds to the formation of the short-scale structure beyond the scattering 
time scale within a window comparable with the interaction range. The second inequality expresses 
the persistence of this short-scale structure within this spatial window up to times much longer than 
the scattering time. 

Remark 4- Note that in the definition of .F/v(t) we compared ^N,t/0- ~ ^n) with the initial 
product state and not with the evolved product state cpf which would have been more natural. 
For the relatively short time scales that we can consider, this difference is irrelevant; the condensate 
essentially does not move. One can directly check that J 9tix\ — X2)\<p® N i~x) — (/9f ljY (x)| 2 dx <C JFtv(O), 
hence the modification does not influence (|1.21j) . To investigate the persistence of the local structure 
up to times of order 1 in N, the definition of ^jv(t) should, of course, contain (pf N instead of the 
initial state ip® N . However, at such large times, the i-th and j'-th particles also interact with other 
particles. Our analysis does not control consecutive multiple collisions, although the validity of the 
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Gross Pitaevskii equation still gives an indirect evidence that the correlation structure is preserved 
even for times of order one. 

Remark 5. In [2] the local structure was identified by the L 2 norm of the mixed derivative 
ViV2[ x I / Af,t/(l — &n)]- For initially factorized states, the integral 

*jv,t(x) 2 



9 e (xi - x 2 ) 



ViV 2 : 



dx 



(1.23) 



'1 - u N {xi - x 2 ) 

is of the order N at time t = 0, and it is expected to be of order one for times t S> iV~ 2 . Proving this 
decay would establish the formation of the local structure in a much stronger norm than the local L 2 
norm used Theorem ll.il Unfortunately, due to the very singular interaction potential Vn, we cannot 
bound (|1.23p effectively (estimating it by the expectation of produces a bound proportional to 
N for all times, see Lemma iB.ip . 

Remark 6. Condition (j!.18[) encodes all regularity and decay that we assume on <p, although it is 
not optimal and we do not aim at finding most general conditions on (p. In particular, (|1.18p implies 
that (p G H 3 (R 3 ), hence <p G C^R 3 ). 

Notation. By C we will mean a constant depending only on the unsealed potential V and the 
initial one-body wave function ip. Constants denoted by are meant to depend also on the 

indices p, q, etc. 



2 Proof of main Theorem 

In this Section we present the main steps of the proof of Theorem 11.11 

Proof of Theorem \l.l[ Let ^jv(i) be the quantity defined in (11.16|) . 

To evaluate Tjy, we introduce a dynamics where particles 1 and 2 are decoupled from the others. 
We define 



h (l,2) 



H 



(3) 
N 



u 



(1,2) 



N 



-At - A 2 + V N ( Xl - x 2 ) 

N 

^(-Ai)+ VN{xi-Xj) 

i=3 3^i<j^N 
N 

{Vn(xi - Xj) + V N (x 2 - Xj)) 

3=3 



H 



Then we have 



N 



iH N t r J$N 



. (1,2) rr(3) 



(1,2) 



Hn — U N - 



t ■ 



with ifa = e-^tip® 2 and $ t = e ~iHft^(N-2)^ TYm ^ 



^jv(t) < / 0t(xi-x 2 ) 



+ I 9e(xi - x 2 ) 

< c(g N (t)+Jc N (t)) 



^jy,t(x) - ^jy,t(x) 
1 - uj n (xi - x 2 ) 

*iv,t(x) 



dx 



N 



1 - u N (xi - x 2 ) 



1=1 



,1 i 



dx 



(2.1) 
(2.2) 

(2.3) 
(2.4) 

(2.5) 



(2.6) 
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where 



G N (t) 
IC N (t) :-- 



9^(xi - x 2 )|*jv,t(x) 



9 e (xi - x 2 ) 



*jv,*(x) 



*jv,t(x)| dx, 
N 



1 - un{xi - x 2 ) 



i=l 



dx 



(2.7) 
(2.8) 



(notice that 0^(cc) < @ 2 £(x) has been used). 

£/jv measures the L 2 -distance, in the spatial region where \x\ — x 2 \ < £, between the evolutions 
^N,t and ^N,t of the initial state ip® N with the dynamics given by Hjy and Hn, respectively. By 
construction, Gn(0) = and at later times <7jv(i) deteriorates as the two vectors ^N,t an d V&jv,* 
separate. We control such behavior in Section [3) the result (Proposition 13. ip is 



Q N (t) < C(N log N)st 2 



(2.9) 



provided that t <C iV~ 2 / 5 and for all times t ^ 0. Notice that, if we had followed the dependence on 
the length scale A characterizing the initial wave function ip (as explained in Remark 1 after Theorem 
II. ip . the bound (|2,9p would have taken the dimensionally correct form GN(t) < C(iV log 7V) 4 / 5 (i/A 2 ) 2 ; 
similar remarks would also apply to all estimates in the sequel which are apparently dimensionally 
incorrect. 



On the other hand, since factorisation = ipt ® &t is preserved in time, see (j2.5j) . then JCn 
turns out to be essentially a two-body quantity, that is, an integral only in variables x% and x 2 , which 
involves the Schrodinger evolution of the initial datum ip® 2 with interaction Vn(x\ — x 2 ). Reduction 
of /Cat to a two-body integral and a remainder is done in Section [H The result (Proposition I4.1| ) is 



IC N (t) ^ Ct / d V dx9 2e (x) 

Jl 3 xl 3 



-it) N t 



1 - u) N (x) 



+ C 



(log N 2 t) 
WT 



+ t 2 +£ 3 m + e 5 



(2.10) 



for all times t > 0. Here we defined 



t) N = -2A x + V N (x), and ip v {x) = ip(r) + x/2)<p(r) - x/2) (2.11) 

(77 and x denote, in other words, the center of mass and, respectively, the relative coordinate of 
particle one and two). 

Last, we consider the two-body integral contained in our bound (|2.10|) to /Cat. It has the natural 
interpretation of a quantity which tracks the dynamical formation of a short-scale structure in the 
evolution of the two-body initial factor state <p® 2 . We study this problem in Section [5j The result 
(Corollary E2]) is 

/ drjdx9z(x) 
Jr 3 xR 3 

for all times t > 0. 

Bounds (I23D, piDjl . and (I2T21) complete estimate for F N . It takes the form 
F N (t) < c[{N\ogN)U 2 + ( lo ^*>V + t 2 + l 3 iV* + l 5 + {l ° g ^ 2 2 t t)6 N 2 £ 5 +£ 2 Nt 2 ^ (2.13) 



e- lt)Nt ip v ){x) 
1 - uj n (x) 



(log iV 2 t) b 



7V 2 £ 3 + Nt z + f 



(2.12) 
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for t > and iV _1 < £ < iV _ 5. Using that, by Lemma Ell F N (0) ~ £/N 2 and restricting to 
£ <C iV" 1 / 2 , f < AT -1 , one has 



^ (t)< C^(0),<^<^+<^ (1 o g ^)« 



(2.14) 
□ 



We compute now the asymptotics of at time t = 0. 



Lemma 2.1. Lei .T-jv fre the quantity defined in 111.16)) . Assume that £ scales with N in such a way 
that N£ — * Co when N —* oo, where Co G (0, +oo]. TTien i/iere exists a constant C x , depending only 
on the cut-off function x defined in $1.15]) and on the potential V, such that 



N 2 

hm —^(o) = c x y\\i 

Nl-+co 



Moreover, if N£ — > oo, 



hm 

M-*oo 



7NP 



■^iv(O) = 4vra 2 ||xlUi(K)l|y|ll4 



(2.15) 



(2.16) 



y9(xi)y9(x 2 ) 
1 - wjv(xi - x 2 ) 

Wjv(£ 



Proof. One has 

AT 2 AT 2 /• 

—^(0) = — / 9 l {x l -x 2 ) 
£ £ Jw 6 

N 2 f f 
= —— I dx 2 \p(x 2 )\ 2 / dx6 e (x 

= ^/ M 3 dx2| ^ 2)|2 l dxx (M) 

TV— >oo ^ II ||4 

by continuity of ip and by dominated convergence, where 



(p(xi)ip(x 2 ) 

2 



1 - Wjv(x) 

1 — uj{x) 



\<p(x 2 + x)\ 



dxidx 2 

2 



.r 



C v : = lim 



o;(x) 



1 — 

The above limit clearly exists if cq is finite. If, instead, N£ — » oo one has 



1 

iV7 



1 — u;(x) 



2 Nl-^oo , 2 n n 
► 4 ™ ||x||ii( 



To prove (I2.19p . one sees that integration for \x\ > 1 gives the leading contribution, since 

x Ni-too 



1 

m 



\x\<l 



Uj(x) 2 _ 1 47r 



1 - u(x) N£ 3 (1 - cp) s 







(2.17) 



(2.18) 



(2.19) 



(2.20) 
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while, when \x\ > 1, since to{x) = a/\x\ 



<?[ dx x(M) 

NiJ\ x>1 (\x\-ay 



Ana 2 I drx( 



2 Jl 



H (2.21) 



4™ 2 Hxllii(M) 

and (|2.16|) is proved. □ 

3 Many-body problem for short times 

In this section we control the growth of the nonnegative quantity Qn defined in f)2.T|) . 

Proposition 3.1. Assume that I <C N~ 2 / 5 . Then, 

G N (t) ^ C(N log N)tt 2 (3.1) 

for all times t ^ 0. 

Proof. For arbitrary i > 21 we have 

&r(t) = / ^(xi-x 2 )|*iv, i (x)-$iv,t(x)| 2 dx 

^ / ^(x 1 -x 2 )|^ 7V ,t(x)-$iv,t(x)| 2 dx =: Q N (t). 

JR3N 



(3.2) 



The parameter £ will be fixed later on (we will choose £ = (iVlog N)~ 2 / 5 ). Let us denote by (•, •) the 
scalar product in L 2 (R. 3N ) and by #12, uj\ 2 , and Vij the operators of multiplication by 9t(x% — x 2 ), 
ojn{x\ — X2), and Vn(x% — Xj) respectively. Then 

= ^N,t ~ *N,t , #12 i*N,t ~ *N,t)) 

= (V N , t - &N,t , [iH N , e 2 l2 ] (^ N>t - *iv,t)> (3.3) 
+ 21m < 9f 2 (^ N>t - #N,t) » (H N - H N )y N>t ) 

Since 

[H N ,6 2 12 \ = [-A 1 -A 2 ,6 2 12 ] = (v r -(V^i 2 )^ 2 + 12 (V r ^ 2 )-V r ) , (3.4) 

r=l,2 

the summand on the r.h.s. of (|3,3p takes the form 



2i Yl (( V rMN,t,(Vr0l2)0 12 5y Nit ) - ( 5V Njt , 12 (V r 12 ) V r (5*^)' 



r=l,2 



(3.5) 
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having set 5^N,t '■= *&N,t — ^N,t- Then 

I Jn® I < 4 S I (( V ^12) • (Vr<5*JV,t) , #12 6* N ,, 



r=l,2 



< 8\/&K*) ( / dx|V^i 2 | 2 |(Vi<5^)(x)| 2 



Q N (t) IIW12IL ( / dx|(Vi<5* A r, 4 )(x)| 2 



2 



Let us now examine the summand in the r.h.s. of (13.30 . One has 



JjSF J (*)| < 2|<6>? 2 (^, t -^),(ff JV -fl JV )* 
= 2|(<5^,^ 2 ^ 1 ' 2 )$ 7Vit ; 

JV 



JV.t/ 



,t 



^ CyJg N (t) (N(y N!t ,ef 2 v? 3 y N>t ) + N 2 (y N!t ,el 2 Vi3VuVN : 

The first term on the r.h.s. is estimated by 

(^ N)t ,6j 2 V? 3 y Ny t) = / 9j{xi - x 2 )V${x l - x 3 )\Mxux 2 )\ 2 \$t{x 3 , ■ ■ ■ ,x N )\ 2 dx 

' - n v 



< UtWto 6 z i (x l -x 2 )N li V\N{x 1 -x 3 ))\<5> t {x 3 ,...,x N )\ z &x 
= CUtWloN^f V 2 (Nx)\<S> t (x 3 ,...,x N )\ 2 dxdx 3 ---dx N 



/R3(iV-l) 



c\m\l>Np\\*t" 2 



(3.6) 



< c^g N (t) r 1 . 

On the last line above we used 

/ dx|(V 1( 5^ t )(x)| 2 <C. (3.7) 

71R3JV 

Indeed, by the symmetry of ^N,t, and because of (jB.ip . 

/ dx|Vi*^ >t | a = (*jv,t,(-Ai)*jv,t) < (^ w ,^0 <C, (3.8) 

J K 3iV iV 

while, due to the factorization (|2.5p . 

dx|Vi*iv lt | 2 = f dxidxaKVie-^V^Cxi.xa)! 2 < <¥>® 2 , f)£'V 2 ) < C. (3.9) 



(3.10) 



< 4|(^,^ 2 (^y lj )^, t )| 

J=3 

(last line following by the permutational symmetry 1 <-» 2), whence, by Schwarz inequality, 

i / N \§ 
|^ 2) (t)| < 4<^ JV ,*,0? 2 ^jv,t> 5 I (^m^nVikVij^^t)) 

V i,fc=3 ' (3.11) 



(3.12) 
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It follows by Proposition IA.1I that, under the assumption (|1.18p . 

ll^li~( K 6, dxidX2) <C(logA0 (3.13) 

and thus that 

(^ N ,t,0 2 12 Vi 3 ^N,t) ^ CN(logN) 2 £ 3 . (3.14) 
On the other hand, the second term on the r.h.s. of (|3,lip is estimated as 
(V N>t , el 2 V X zVu ®N,t > = 

a2/„. „ \ rr_/„ „. \t/:_/_ ... \ L;. f„ „ \|2 Ij, /„ \|2 



^(x! - x 2 )V N (x 1 - x 3 )V N (xi - x 4 )\x(jt{xi,X2)\ \® t {x3, - ■ ■ ,x N )\ dx 
< c \\tpt\\lo^ 3 V N (xi-x 3 )V N (x 1 -x 4: )\®t(x3,---,x N )\ 2 dx 1 dx 3 ---dx N 

7R3(W-1) 

^ 3 ||Vjv||3 / dxs-'-dxAr 1^1 - A 3 $i(x 3 ,...,a;jv)| 2 / do* yjyrfci - 

2 J K 3(JV-2) J R 3 



< c 



|2 ?3| 

loo ' 



2 



(3) 



< c Uttc WhWh-^r (* t , ( i + ) «T> 



TV \ V N -2J /i2(K3(JV-2)) 



In the last inequality we used again bound (|3.13p and the asymptotics (jB.ip . Thus, (|3.1ip reads 

\j£\ t )\ ^C^G N (t) Nil log N. (3.16) 

Altogether, ([331), flESH, and ([3T6]) give 

|^v(*)| < C^g N {t) (r 1 +NlhogN) . (3.17) 
Letting £= (NlogN)~ 2 / 5 , we get 



2 



|&r(<)| < Cy/G N (t) (iV log AT) 5 (3.18) 
which implies (|3.1|) by Gronwall Lemma, because <7jv(0) = 0. □ 

4 Reduction to the two-body problem 

The goal of this section is to reduce the study of the quantity ICn, defined in (|2.8p . to the analysis 
of a two-body term (which will then be controlled in Section [5]). We will use, in this section, the 
coordinates (77, x) defined by 

77 = (x\ + X2j/2 (center of mass coordinate) ^ 
x = X2 — x\ (relative coordinates). 

In these coordinates, the two-body Hamiltonian t)\r' introduced in (|2.ip takes the form 

h^' 2) = -Ar,/2 + \) N , with t) N = -2A x + V N (x). (4.2) 
Note, also, that the two-body initial data ip® 2 is given by 

i/j(ri,x) = ipr,(x) = ip(ji + |J <p(rj- |) . (4.3) 
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Proposition 4.1. Suppose that the assumptions of Theorem \l.l\ are satisfied. Let fj^r 
and il>rj(x) be defined as in Then, if K,N{t) is defined as in \2. 6\) . we have 



-2A+V N {x) 



K N (t) ^ Ct 



d.77 dx 2 £ (x) 



V,: 



1 - uj n {x) 



+ C 



( logNH)6 ^ + t 2 + £ 3 Nt + £ 5 ) (4.4) 



NH 



for all times t > 0. 

Proof. We divide the proof in four steps. 
Step 1. We have 

IC N (t) < 2 / 9 e (xi-x 2 ) 

J«?xM. 3 



ifj t (xi,x 2 ) 
1 -wat(xi -x 2 ) 



ip(x 1 )ip(x 2 ) 



dxidx 2 + CrWt 



where tpt(xi,x 2 ) = exp(— it\)^' 2 ^)(p® 2 is defined in (12. 5h . 
To prove (|4.5p , we use that ^ at * = <8> &t and we split 



/Civ(t) 



9 e (x 1 - x 2 ) 



N 



1 - uj n (xi - x 2 ) 



i=l 



dx 



^ 2 



9 e (xi - x 2 ) 



A(xi,x 2 ) 



1 - ui N (t)(xi - x 2 ) 



ip(x 1 )ip(x 2 ) 



dxidx 2 



+ 2 



e t ( Xl -x 2 ) (^ 2 ®^)(x)-^(x) 



dx 



where 



dx 



? e ( Xl - x 2 )|</?(xi)| 2 |(/?(x 2 )|' 



i(N-2) 



L 2 (R3(iV-2)) 



s: C£ 3 

Eq. (|4.5p follows now from 
To show (14.81). observe that 



|$t - $t=o| 



L 2( R 3(iV-2)) • 



d_ 
di 1 



— - $ || 2 



< 2|<d> , J f/}, 3) <I> t >| < 2<cD ,^ 3) $ ) < CN. 



(4.5) 



(4.6) 



dxidx 2 

(4.7) 



(4.8) 



(4.9) 



The last inequality follows from (jB.ll) with N replaced by N — 2 (recall the definition of Ilffi from 



Step 2. We have 
0<(xi - x 2 ) 



ipt(xi,x 2 ) 



1 - uj n (xi - x 2 ) 
<C£ 2 I dj]dx9u{x 



<p(xi)(p(x 2 ) 



dxidx 2 



1 - LO N (x) 



+ 



c(nV(t)+K%\t)) 



(4.10) 
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with 



11$ (t) := I drjj dx9 e (x 



dy9 e (y)- 



and 



11$ {t) := I dril dx9 e (x) 



Here we use the notation 

To prove (|4.10p . we observe that 
dxidx2 Oi(x\ — X2) 
drj dx 9e(x) 

<R 3 

drj dx 9i(x 



1 - u) N (y) 

2 iP(ri,x)- dy9 e {y) r - — '- 

Jr 3 1 - u N (y) 

0t(v) 



t{y) 



ipt(xi,x 2 ) 
1 - u N (x 1 - x 2 ) 



(4.11) 

(4.12) 
(4.13) 



ip(xi)(p(x 2 ) 



e i-%Lt ( e -^^)( x ) 

1 - u N {x) 



ip(ri,x) 



«S 2 



e % 2 * (e % 



+ 2 1 drjdx9 e (x 
'l 3 xt 3 



^ 2 



drj dx 9i(x 
+ 2 dr]dx9 e ( 

./R 3 xR 3 



1 - uj n (x) 

ip(r],x) - / dy9 £ (y) 
Jr 3 

(e- il > Nt tp v )(x) 



.3 l-u N (y) 
e i^t ( e -i^*^)(y) 2 



(4.14) 



1 - wat(x) 



dy ^(y) 



1 -u N (y) 
(e-^^iy) 2 



t ip(i],x) 



e 



dy 9 e (y) 



1 - wjv(y) 



1 -u N (y) 



by unitarity of e lA,)<//2 . The second term can be clearly bounded by the sum of 11$ (t) and K${t). 
For every fixed ?] 6 I 3 , the first term on the r.h.s. of the last equation can be estimated using the 
Poincare inequality as 



/ dx9i(a 



(e~*^Vij)(a5) 



dy 9 e (y)- 



il,Nt Tpn)(y) 



1-oj n (x) 7 M 3 a " Lyaj l-uj N (y) 
(e-^ Nt rp v )(x) 



I dx 

'supp e t 

< CI 2 f dx 

J SUpp 0£ 



, d^fo) ( e " to '-wfa> 

(e- if > Nt i/> v )(x) 2 



V x - 



(4.15) 



S'iep 3. Suppose that 11$ (t) is defined as in (I4.1ip . Then 



V 



1 - uj n (x) 

(e-^ Nt ip v )(x) 
1 - uj n (x) 



(4.16) 
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To show (|4.16p . we note that, from (|4.11|) . 

n%\t) < ce 3 J^ d v || e -*^* w ^ H^^s^) • (4.i7) 

Let J7jv be the wave operator associated with the Hamiltonian f}jv, defined as the strong limit 

Q N = s - lim e it,Nt e 2iAt . 

t— >oo 

Then, by the intertwining property (jD.ip and Yajima's bound (|D.2p . 

||e-^W^t < C||e 2i * A ^^^||L- (4.18) 
In Proposition 15.31 we prove that 

||e 2i * A OW^||L < ^3/, V S E(3,+oo] (4.19) 

where c s ~ (s — 3)~ 6 as s — ► 3 + and where we defined the norm 

IWI = ll^llvya,! + ll^llwrs,^. (4.20) 



Optimizing in s > 3, we have 



and thus, since \\ip v \\ < C(rj)~ a for some a > 3 (by the definition (|4.3p . and the condition (|1.18p ). 
we obtain (14.161) . 

Step 4- Assume that 1Z N (t) is defined as in (|4.12p . Then 

T^n (*) < C (* 2 + ^ + • (4-22) 

First we rewrite 

dygM ^^'T^^ = / dy^( y )(e- 2 ^^)(y) (4-23) 
i - ^Af(y) Jr3 

by means of the operator 

£jv := -A + 2 n VcJjV -V. (4.24) 
1 - Wat 

In fact, since P)at(1 — wat) = (— 2A + Vat)(1 — ^n) = 0, one has 

1 



1 - Wat 

whence 



hAr(l -waO0 = 2£ N <p V<p£H 2 (R d ) (4.25) 



e - if )N* (1-^)0 = (l-u N )e- 2i£Nt (j) V0G L 2 (M 3 ). (4.26) 
It is worth noticing that 

(</>, £^). = (£ N <f>, V). = (V0, VV). (4.27) 
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where 



(/, g)» = j da; (1 - lo n (x)) 2 f(x) g(x) . 



It follows that the operator Cn is self-adjoint on the weighted Hilbert space L 2 (M 3 , (1 — uj]y(x)) 2 dx) 
{Cn is the Laplacian on the weighted space). It is also important to note that, because of the 
properties of ujn, the norm || • |. defined by the weighted product (•,•)» is comparable with the 
standard L 2 -norm, in the sense that c||</>||2 < \\4>\\» ^ ll^lb, with an appropriate constant c > 0. 

From ([4TT2]) . we find 



+ cf d v [ dxe e (x)\^ v (x)-^(o)f 

JR 3 JR 3 



+ C dr, dx9 e (x) ^(0)- / dye e (y)ip v (y) 



+ C dr?/ dx6 e (x) / dy9 e (y)(e 



-2itC r 



i)^(y) 



(4.28) 



To estimate the first summand on the r.h.s. of (|4.28p we use that 



— I 

dt 



e * 2 ip(v, x ) ~ VK 7 ?' x ) 



\L 2 (R 3 ,dr]) 



< 2||V^(77,^)||i2 (R 3 )d ^ 



^ C 



uniformly in x £ M 3 (the last inequality follows from (|A.18|) ). whence 



(4.29) 



dx0t{x) / dr/ e * 2 t ij}{rf, x) — tp(r,, x)\ < Crt. 



R 3 



(4.30) 



To control the second summand on the r.h.s. of (|4.28j) . we observe that, 



ifj n (x) - tp v (0) 



f 1 d 

/ ds—ipnisx) 
Jo ds 



^ \x\ ||V x ^n||oo < C- 



(vY 



(the last inequality follows from (|A.18|) ). and thus 



(4.31) 



dry/ dx0 e (x)\^r,(x)-iP v (O)\' ^ C 



dr) 



(v) 



2a 



dxei(x)\x\ 2 < ce 5 . 



(4.32) 



The third summand on the r.h.s. of (|4.28p can be bounded by 



dr, dx8 £ (x) ^(0)- / dy9 e (y)^ v (y) 



dr, dx0 e (x) / dy6z{y){i> v {y)-^)) 



< C 



dr. 



(V) 



2a 



dx6 e (x)[ / dy9 e (y)\y\ ^ Ct 



(4.33) 



where (14.311) has been used. 
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Finally, to estimate the fourth summand on the r.h.s. of f|4.28[) . we observe that, for fixed rj E M 3 , 



dye e (y){(e- 2ltCN -l)^ v )(y) 
Expanding 



-2itCr 



-2itC? 



1)^|| 2 (4.34) 



-2it£jv 



l)i> v = -2i dse- 2isCN £ N ip v 
Jo 



we obtain 



dyO t (y)((t 



i)Vy)(y) 



2 ct r, „ 

-Wo" 5116 



<^ j\s\\e- 2is ^C N ^\\l (4.35) 



<r Ct \\r / iia 



To bound the r.h.s. of the last equation we observe that 

||^V"0t?||. ^ 11^^11+2 Vuin ■ Vi/'j? 



< C(7?)- Q + y dx|V^(x)| 2 |V^(x)| 2 



(4.36) 



for a > 3. Here we used the definition (ET24T) of £ N and Eq. (OS]) . Plugging (|4T34) and (|4T36l) into 
the fourth summand on the r.h.s. of (|4.28p . we find 



drW dx6 e (x) / dy^(y)((e~ 2 ^ JV -l)^)(?/) 



< Ct 2 



The results of Steps 1-4 complete the proof of (|4.4p . 



(4.37) 



□ 



5 Dynamical formation of correlations among two particles 

In this section we estimate the quantity 



d-q dx 02e(%) 



V,: 



1 - LO N (x) 



(5.1) 



which arises in the bound (|4.4p . 

To control the integral (|5. 1H . we are going to make use of the following proposition, which is 
stated in terms of macroscopic coordinates. 

Proposition 5.1. Suppose that V is a non-negative, smooth, spherically symmetric, and compactly 
supported potential. Let f) = — 2A + V and denote by 1 — u the solution to the zero energy scattering 
equation 

b(i-w) = o 
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with boundary condition uj(X) — > as \X\ — > oo. Moreover, let 9^ be defined as in \1. 1$ for some 
L > 1. Consider ip £ W 3 ' 1 ^ 3 ) n W 3 '°°(M 3 ) ; wtfi ||^|| = 1. T/ien, /or A > 1, de/me ^ A as 



^a(X) := VWA) 



TOf/i A>I. Define 



Then 



Fk,l(T) := / 



V>a)PO 



1 - 



dX. 



for all T > and A>I sufficiently large. Here we used the notation 



W3,l + ||l/)||w3,oa 



(5.2) 

(5.3) 
(5.4) 

(5.5) 



Remark 1. It is simple to check -Fa,l(0) ~ 1+L 3 /A 2 . Thus, Prop osition 1 5 . 1 1 st ates that F\^{T) <C 
^a,l(0), if L 3 <C A 2 , and times T such that 



1«T«A. 



(5.6) 



This fact can be interpreted as a sign for the formation of the local structure 1 — u> in e l ^ T tp\. 

Remark 2. Taking formally A = oo, (|5.4p describes the relaxation of a constant initial data 
towards the solution 1 — uj of the zero-energy scattering equation. 

Remark 3. For A < oo the function ip\ can be thought of as cutting off if)^ = 1 at distances 
A. Since the energy hence the velocity of ip\ is of order 1 in A, a time of order A is necessary for 
the effects of the cut-off to reach the window of size L: this explains why we can only prove that 
Fa,l(T) <C F a , l (Q) for times T«A. 

Applying Proposition 15.11 we immediately obtain the following bound for the integral (|5.ip . 
Corollary 5.2. Under the assumptions, and using the notation introduced in Proposition \4 ■ 1\ 



dn dx #2f(x) 



V, 



(« 



-it) N t 



Ipr,) (x) 



1 - uj n {x) 



^ C 



(log NH 
NH 



N 2 e 3 + m z + . 



(5.7) 



for all times t > 0. 



Proof. Changing coordinates to X = Nx, we obtain, from Proposition 15.11 (with L = 2N£ and 
A = N), 



dx 2 e(x) 



-it} N t 



<C\ 



,f( {l0g ^ f NH^Ne + A. (5i 



1 - u N (x) 

Since, by (1A.18|1 . IV^I < C (rj)~ a for some a > 3, (15.70 follows from (|5.8p . integrating over rj. 



□ 



Remark. By Corollary 15.21 the integral (|5.ip can be shown to be smaller than its value at 
time t = 0, for all times t in the interval N~ 2 <C t <C ./V -1 . Note that in the many body setting of 
Theorem ll.il on the other hand, we can only prove that ^ r /v(t)/^ 7 /v(0) <C 1 up to times t <C A^*- 2- !") 
(if £ ~ A^ 1 ); this is due to the lack of control of the many body effects for larger times. 
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Next, we prove Proposition 15.11 As we will see, the two main tools in the proof are Yajima's 
bounds on the wave operator f2 associated with the one-particle Hamiltonian —A + and a new 
generalized dispersive estimate for initial data which are slowly decreasing at infinity but have some 
regularity. The dispersive estimate is presented in Section [6J The definition and the most important 
properties of the wave operator are collected in Appendix [Dl 



Proof of Proposition \5.1\ We start by splitting 

F A>L (T) ^ 2*£i(T) + 2F| 2 [(T) 

where 



F ( aI(T) := [ 9 L {X) 
F^JT) := f 9 L (X) 



-if)T 



V 



V 



to 



V>a)PO 



dX . 



1-uj(X) 
-^ t (1-^a)(A) 



l-u{X) 



dX. 



(5.9) 

(5.10) 
(5.11) 



Let us first estimate F^ L . Let Q, be the wave operator associated with —Ax + \V{X) (that is, 
our ^()) as defined in Proposition ID. 1[ Then, by the intertwining relation (|D. 1|) . 



F£l(T) = [ dX9 L {X) 



V 



(Fie 2iTA WluiP a )(X) 



1-lo(X) 



(5.12) 



< c \\Vuj\\ 2 2 \\ne 2iTA n*uj^ A \\ + l 3 || vne 2iTA n*uji/jt 



Here we have also used that 1 — oj{X) ^ const > (see Lemma IC.lj) . By Yajima's bound (]D.2jl . we 
have 

(5.13) 
(5.14) 



rJt lloo II Tli Moo 



and 



VOe 2iTA Q*^ A \\I^C( \\ e 2tTA fi*u^ A £ + || Ve 2,rA 0*o^A £ 



Therefore, using (|C.3p to bound ||Vo>| 



, 2«TA r>* 



(5.15) 



L 3 



Vs G (3, oo] 



(5.16) 



From Proposition 15.31 below, we get 

F { A %(T) < c s Uf 

where c s ~ (s — 3)~ 6 as s — > 3 + . 

Let us now estimate the term F\ L defined in (|5.1ip . We rewrite it conveniently by means of the 
operator 

C := -A + 2 V 

1 — UJ 



already introduced in (|4.24j) (in microscopic variables). Analogously to (I4.26h . we have 

e- it)T (1 - u)<f> = (1 - u)e- 2iCT ^ , V 4> G L 2 (R 3 ) 
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and 

{<j),C^). = {£<f>,ip}. = (V<£,VV). 

where bulleted scalar products here are in the Hilbert space L 2 (M 3 , (1 — uj{X)) 2 dX\ . On this space 
C acts then selfadjoint operator. So 



= I dX9 L (X)\Ve~ 2lCT i; A (X)\ 2 
Jk 3 

^ [ dX6 L (X)\ V^ A (X)\ 2 + [ dX|V(e 



The first summand in the r.h.s. of (|5.17p is bounded as 



dxe L (x)\Vip A (x)\ 2 ^ c^\\vnlo ^ 



-2iCT 



L 3 



A 2 ' 



The second summand in the r.h.s. of (|5.17p is bounded by 



dX\ V(e 



-2iCT 



1)?Pa(X) 



T 

V/ dSe 



-2iCS 



CM*) 



^ C dX(l -Lo(X)f 



V / dSe 



-2iCS 



^ C 



/ dSVe- 2tCS C^ A 
Jo 



< CT A sup || Ve~ 2iCS £4j A \\' 

Se[o,r] 

< CT 2 ||V£^a|^ 



because, by (|4.27|) . 



Ve- 2i£5 £VA ||: 



(Ve- 2i£5 £VA,Ve- 2i£5 £^A>. 
{e- 2iCS C^ A ,Ce- 2iCS C^ A ). 

{Cip A ,C 2 iPa)» 
(V£Va,V£^a). 

|| V£^ A ||! • 



In turn, || VCtp A \\ 2 is estimated as 



V£Va||. ^ 2|| VVa||. + 2 



^ 2||V 3 Va||^ + 2 



, Vw 
V ( VVa 

1 — U} 

1 — w 



where 



V 3 ^a| 



v 3 Vil 

A 3 



A 3 



and 



) |[ < C ( || |V 2 w| |Wa| IIJ + II |V^| 2 |VVa| K + II |Vw| |VVa| H2 ) 



A 2 



(5.17) 



(5.18) 



(5.19) 



(5.20) 



(5.21) 



(5.22) 



(5.23) 
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Thus, 



/ 



dX\V(e- 2 * CT -l)MX)\ 2 < C^IWI 2 . (5.24) 



Altogether, we find 

T 

j?( 2 \(t\ <r n III ,/Jl 2 

A 2 

Finally, plugging estimates (|5.16j) for F^l and (|5,25p for i^? 2 ]- into (|5.9p . one gets 



F| 2 i(T) < C|W 2 ^. (5.25) 



Fa,l(T) ^ C HI VIII 2 (c*^ + r2 ^" 2 L3 ) Vs G (3, oo] . (5.26) 
Optimizing in s (since c s ~ (s — 3) -6 as s — ► 3 + ), we find 

Fa,l(T) < C|||V||| 2 (i^L 3 + (5.27) 

and (|5.4p is proved. □ 



The following Proposition, which played an important role in the proof of Proposition 15.11 and 
was also used in Step 3 of the proof of Proposition 14.11 provides an estimate for the decay (in L°° ) of 
the evolution generated by the Hamiltonian t) = — 2A + V(x) on initial data decaying only as |x| _1 . 
It is based on new dispersive bounds for the free evolution which are presented in Section [6J 



Proposition 5.3. Under the same assumptions, and using the same notation as in Proposition \5.1 
we have, for every T € M, 



\\e- iTA n*uip A \\ wl ^ ^ c s ^l VsG(3,+oo], (5.28) 

uniformly in A. Here c s ~ (s — 3)~ 3 as s — ► 3 + . 

Remark. For our purposes, the bound (|5.28p is better than the standard L 1 — ► L°° dispersive 
estimate, because it is uniform in A and we eventually want to let A —* oo (on the other hand, since 
to only decays as |x| _1 , the standard L 1 — » L°° dispersive bound would diverge like A 2 ). 

Proof. Let be the wave operator associated with —Ax + \V(X). We will use the dispersive 
estimate (|6.4p for s G (3,oo], which gives 

|| e~ lTA WluiPa \\ 2 ^ ( || n*u^ A \\ 2 s + \\ Vfl^A fju + II V 2 ^Va f_3^ ) (5.29) 

_/ / \ s+3 2s+3 / 

and 



Ve- iTA ^^ A ||L < -^(HV^VaII" +||V 2 ^^ A ||k + ||V 3 ^VL|| 2 _3^), (5.30) 



respectively. For convenience we set 



3s / 3 

u:= 7T3H2' 3 



3s (a 3 

and z := G 1, - 

2s + 3 V ' 2 



(5.31) 
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We shall make use of the following bounds, which are a consequence of the properties of uj (see 
Lemma IC.ip . and are proven separately in Lemma 15.41 



||V m MA)|| P < Cp, 

|v m (y^ A )|| p ^ c 



Vp G [1, oo] , Vm G {0, 1, 2, 3} s.t. p(m + 1) > 3 (5.32) 
Vp G [1, oo] , Vm G {0, 1} . (5.33) 



One sees that Q*uj4>a £ £ S (M 3 ) Vs > 3 because, by Yajima's bound (jD.2[) and by estimate (|5.32p . 

||n*w^ A || a < c s |Ma||s < c,|^|. (5.34) 

Here c s ~ (s — 3) -1 / 3 for s — > 3 + . 

To treat || V£1*ujiJja\\v we use the property 

||V/|| p < cp || V=A/|| p VpG(l,co) (5.35) 

valid for every / G -L 2 such that \/— A/ G L 2 (see Lemma l5.5f) . This is the case for / = O*w0a, 
due to Yajima's bound, although the corresponding norms are not bounded uniformly in A. In our 
application we need p G (f,3], so we are never at the borderline case and c p remains uniformly 
bounded. Thus, we get 



V£l*u>ipA || ^ || v 7 — A £l*u}ipA 

^ Cv 



c v fJ*yf) ojipA 



rj ojipA \\ v < Cv 



dk 



5^ C)! 



dk 



o Vk 



t 



-A + k 



o Vk 



1 



k + t) 



(5.36) 



where we have used intertwining relation (|D.1|) . Yajima's bound (|D.2|) . the fact that Ij generates a 
positivity-preserving semigroup, and 



1 



dk 1 



We notice that estimates ()5.32|) and ()5.33j) give 

|| fj cxjVaIIo < || — A (oj^a)IU + II^^aIU < c a 
Thus, for any b G [1,3) Pi [l,f), Young's inequality gives 



a G (1, oo] 



(5.37) 



(5.38) 



1 



-A + k 



\t)uipA\ 



< c 



1 

47T 
v.b 



dY- 



-Vk\.\ 



-Vk \X-Y\ 



\X-Y\ 



V 71 



v b ' 



(5.39) 



k2 26 1| fjw^Alln+l-lv 

* D fa ' 



We use this bound to control the r.h.s. of (I5.36p . applying it with two different values of b in the 
region \k\ < 1 and \k\ > 1. For v G (1,3), we find 



(5.40) 



\7n*UJ1p A \\ v < Cv HI frtWA.ll 3v2 + || t)UljjA\\ v 
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Here c v ~ (v — |) 2 when v —* | + For t> = 3 we find, on the other hand, 

|| VO^VaIL < — — fll &<^aII_s_ + II fjo^Alla) (5-41) 

for any // 6 (0, 1). Hence, by f|5.38|) . 

||Vn*u^ A ||^_ < c s U\\ (5.42) 

3+s 

uniformly in A and Vs £ (3, oo], and c s ~ (s — 3)~ 3 as s — > 3 + . 

To treat || V 2c I*loi/j a \\ z we use the C alder on- Zygmund inequality 

||V 2 /|| P < c p \\Af\\ p Vp6(l,oc) (5.43) 

valid for any compactly supported / in W 2,p (W i ) (see [3], Theorem 9.9). The behaviour of the 
constant c p as p —* 1 + can be computed from the constant in the Marcinkiewicz interpolation 
theorem (see [3], Theorem 9.8) and one has c p ~ (p — l) -1 . Although supp(f2*u;?/>A) is not compact, 
we can apply (|5.43p to the compactly supported function % J1*cj^> A ) where we have introduced the 
cut-off function 



at the scale A, with x defined as in (|1.15p . Then one has 

\\v 2 n*u^ A \\ z ^ ||v 2 x a o*^a|L + || v 2 (i - xjn^ A |L 

< c z \\Ax A n*uip A \\ z + \\V 2 (l-x A W^ A \\ z 

< c z \\An*ujip A \\ z + c z \\A{l-x A W^k\\ z + || V 2 (1-x a )^*^a| 

< c|| An*w^ A || a + c z \\v 2 {i-x A W^ A \\ z 

and 



(5.45) 



lim ||V 2 (1-xJ^Va|| =0 (5.46) 

A— >+oo A z 



because 



v\\ - x A )Si"u,p A \\ z < || (v' x Jsr^p A \\ z + 2 1| (VxJ ■ vs!*<Wa||, + 1| (1 - xJvW^aII, 

< c iKrvfcfe + c \m^k + J (1 _ xJv2si .„ fa | L (5 . 47) 

> u . 



Last summand above, in particular, vanishes as A — * +oo by dominated convergence: in fact we 
have (1 — x a )V 2 0*u;t/a — * pointwise and is in L Z (M 3 ) uniformly in A because 

|| (1 -X a )V 2 0*cjVa|| 2 || V 2 n*Loifj A \\ z s$ || 0*^11^2,, < c z \\utl; A \\ w2tZ <oo. (5.48) 

Thus, 

|| V 2 n*u-(p A \\ z < c z || A0*wVa|| 2 V«6(1,+oo). (5.49) 

Using (|5.49p . the intertwining relation (jD.ip . and Yajima's bound ()D.2p . one has (with z = gjxg, 
s>3) 

|| v 2 ^*wVa|| 2 ^ Cz || a^*wVa|| 2 < c z || n*\)ui) A \\ z 

^ c z (\\A{u^ A )\\ z + \\Vu^ A \\ z ) (5.50) 
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uniformly in A and Vs E (3,oo], where (|5.32|) and (|5.33|) have been used in the last line. Following 
the blow-up of the various constants, we see that in (|5.50p c s ~ (s — 3)~ 2 as s — > 3 + . 

From (|5.34|) . (|5.42|) . and (|5.5U|) . the dispersive estimate (|5.29|) takes the form 

II 2 



||e- iTA fnW, A ||^ Vse(3,oc] (5.51) 

with c s ~ (s — 3)~ 6 as s ^ 3 + . 

Next we treat the r.h.s. of (|5.30[) in analogy to what we have done so far for the r.h.s. of (|5.29p . 
The first term on the r.h.s. of (|5.30p is bounded as 

||VO*^ A || s < C || 0*^11^ ^ c s \\co^ A \\ wlt3 ^ clVll (5.52) 

for any s > 3 and uniformly in A: the second inequality above follows from Yajima's bound (ID.2P 
and the last one follows from estimates (15.321) and (|5.33p . The second summand on the r.h.s. of 
(|5.30p is estimated analogously to (|5,50p and gives 

\\V 2 n*u;ip A \\ v ^ CsU\\ (5.53) 

uniformly in A. Finally, the third term on the r.h.s. of f|5.30[) can be bounded, similarly to (|5.49p . 
by 

3 3 

||V 3 ^Va|L < J2\\ v2d i n * u M z < c^||A^*u^ A || z . (5.54) 

i=l i=l 

Hence, by means of the intertwining relation (jD.ip . of Yajima's bound (|D.2p . and of the estimates 
(^32|) and (1Q3D . we find 



A<9;0*u;Va||, = c z || did*l)ujil> A \\ z 

(5.55) 



< c z \\t)uiJj A \\ wliZ 



< c z ( || A Ma)|| 2 + || Vuip A \\ z + || V 3 Ma)|| 2 + || V(Vuip A )\\ z ) 



uniformly in A. Thus, from (|5.52p . (|5.53p . and (|5.55p . the dispersive estimate (|5.30p takes the form 

|| Ve~ iTA n*UTp A \\l < c s M^ V S G(3,oo]. (5.56) 

Following the blow-up of the various constants, we see that in (|5,56p c s ~ (s — 3)~ 4 as s — > 3 + . □ 
In the following two lemmas, we prove estimates which were used in the proof of Proposition 15.31 



Lemma 5.4. Let V , to, ip, and ip A be as in the hypothesis of Proposition HOI Then, for any m S N 
and any p £ [1, oo] such that p(m + 1) > 3, there exists a constant c Pjm such that 

||V">Va)|| p ^ c Pym U\\. (5.57) 

— - 

Here c V)m blows up as (p(m + 1) — 3) p as p — > ( m '^ 1 ) + . Moreover, for any m = 0, 1 and any 
p € [1, oo], one has 

||V m (^A)|| p < CUi (5-58) 
for some constant C depending on V. 
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Proof. We recall from Appendix [C] that oj satisfies the bounds 

||V n cj|| f , < oo Mq G [1, oo] , Vn G N s.t. q(n + 1) > 3 

with 



|V n u;|| g ~ (g(n + l)-3) « as g -> (^)+ . 



Moreover, ip\ satisfies the scaling 



|V>a|| 9 = A"^^ 



(5.59) 
(5.60) 

(5.61) 



Pick m£fJ and p G [1, oo] such that p(m + 1) > 3: then, keeping into account (|5.59j) and (|5.6ip . one 
obtains by Holder inequality 



|V m MA)|| p < Cp,mf2\\ V ^h 



(ro+1) 



(5.62) 



i/=0 



for some constants c Pim blowing up as (p(m + 1) — 3) p when p — ► (^xy) + . (When i/ = it is 



understood that 



p(m+l) 



oo). By interpolation, 



|V m (u;V>A)||p < C p , m y^ ||V^|| g(m+l) ^ Cp im ( 11-011^3,1 + ||^|| W 3,oo) 



y=0 



that is, we obtain (|5.57p . Eq. (|5.58p can be proven similarly. 



(5.63) 



□ 



Lemma 5.5. Let f E L 2 (M 3 ) such that s/^Af G L 2 (R 3 ). Then 



V/ll < Cp||v/=A/ 



(5.64) 



/or any p G (1, +oo). 



Proof. The boundedness of the singular integral operator (see, e.g., [8], Chapter III, Theorem 4, and 
[S], Chapter 13, Theorem 5.1) implies that 



< C||5|| P . 



(5.65) 



Let <? := V - A /. Then g G L 2 (M n ). On L 2 (R n ) the operator y— A has a trivial kernel, therefore 

1 



V^Af = f. 



Then 



V/ll < c P ||v^A/ 



(5.66) 

(5.67) 
□ 
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6 Dispersive estimate for regular, slowly decaying initial data 

The standard dispersive estimate 

h iAt f\\ q < -^\\f\\s lo'Li (6.i) 



3(|-|) IUIIS g = ^ I €[2,+oo] 



for the free Schrodinger evolution is not suited for functions that decay slowly at infinity. In this 
section we prove a dispersive estimate which holds for /6L S for any s G [|, oo], if additionally some 
L p bound is known on the derivatives of /. 

Proposition 6.1. Let s G [|,oo], q G [max{s, 3}, oo], and r G [1, gq^]. £ei / G L S (1R 3 ) such that 



V/GL 



3s 

s+3 | 



p3\ 



V 2 / G L r (R 3 ). 

Then e iAt f G L"(IR 3 ) and 



(6.2) 



II IL < ~ 37F~TT f 11/11. + l|V/|U ) + , C ||V 2 /||r (6.3) 

for some constant C which is independent of s,q,r. 

Remark. We use this estimate in the proof of Proposition 15.31 and of Proposition I A. 1\ with q = oo, 
s G [3,oo] and r = gr§j. In this case, (I6.3P reads 

II eiA V IL < "x( H/IU + HV/II A. + l|V 2 /|| 33 ) . (6.4) 

£ 2s ^ s-f-o 2s-|-3 / 



Proof of Proposition ^. 1[ It is enough to prove (|6.3|) for / G C^°(R 3 ); then the estimate can be 
extended by a density argument. 



For q > 1, we have 



(4vrt)3/2Vy R3 



dx 



z | a; — y \ 

&ye-^u—f{y) 



i\ 1/9 



(6.5) 



We split the above integral for small and large values of \x — y\: in the latter regime integration by 
parts will provide the necessary decay at infinity. We introduce R > and we define the smooth 
cutoff function 

Or{x) := X (M) (6.6) 
with x defined in (|1.15j) . The following scaling properties of Or will be needed: 



9r 




3 

— -<4m,p R p 


1 ' \ m 


p 


I -Or 
1 ' \ m 


p 


3 

= B m , p R v 


V9 R 
1 ' \ m 


V 


3 

— C m ,p R p 



if mp < 3 
if mp > 3 
Vp ^ 1 , Vm G 



(6.7) 
(6.8) 
(6.9) 
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Inserting the cut-off in (|6.5p . we find 



JAt 



c 



q ^ £3/2 



f\\a< 



dx 



+ 



C 



dx 



dy6 R (x-y)e « /(y) 



dy (1 - #r(x - y)) e~u— f(y) 



t 3 / 2 
= (I) + V) 

where first summand in the r.h.s. is immediately estimated by Young's inequality as 



(I) < ^||/IUI|flH|| (1+ l_l)-a < c q , s f3/2 
To estimate the summand (I'), we use 



1 < s < q. 



i i2 I I o 

e^* - = 2rf- • V„ e^~~ 



2 v ^ 



to write 



(I') 



C 



dx 



\x - y\ 



dy (1 - 9 R (x - y)) f(y) -A-^L • V„ e^" 



|x - y| 



2 V S/ 



Then, integration by parts and since 



.-, x ~~ y 

V . _ 

|x — y| 2 



|x - y| 2 ' 



we bound (J') as 



C 



+ 



+ 



dx 



i\x—y\ 2 x — y 

dy /(y) e - — V6 R (x - y) ■ ? 

\x - y\< 



C 

c 



dx 



dx 



dy (1 - 6 R (x - y)) 



f(y) i\*-v\ 



\x - y| 2 



e « 



t l/2 

EE (//) + (///) + (IV) . 



dy (l - 6 R (x - y)) e «' 



i\x-y\ x — y 



\x - y| 5 



V/(y) 



The term (/I) is estimated by Young's inequality and (16. 9j) : 



9i s 
W 2 " 



v q s ' 



t 1 / 2 



1 ^ s ^ 



To control the summands (III), (IV), we integrate by parts once more. From (I6.12|) . 



V„ 



x - y 



\x — y\ 2 \x — y| 4 ' 
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we find 



dy{l-9 R {x-y))-^-e^ 
\x — y\ z 



-2it\ [ d^l-^x-i/jj/^e^ 
+ / dy (1 - 9 R (x - y)) 



X _ yl ,-Vf(y) (6.18) 



dy (V y 9 R (x - y)) ■ ^ e^-f(y) 



\x - y[ 



Therefore 

(///) Ct 

+ Ctt 



^ Cq, s t2 



dx 



dy (1 - 9 R (x - y)) 



f(y) i\*-v\ 



\x - y|- 



e 4t 



da; 

dx 
1 



dy (1 - 9 R (x - y)) -^—J^ . V/(y) 



x — y 

dye « f{y)v9 R {x-y) 



\x - y\- 



+ IIV/IU 

3+s 



(6.19) 



^■3^9 s> 



V9 R 




■ 3 


(l+i-i)- 1 ) 

v q s * 



-*q,s 



+ V/ \\_3s_ 

3+s 



where Young's inequality and (|6.8p and (|6.9p have been used. 



< a < q 



Analogously, to handle the term (IV), we performs a second integration by parts. We observe 



that 



dy (1 - 9 R (x - y)) e « 



ik-yr x — y 



\x - y\< 



V/(y) 



2ii 



dy (1 - ^(* - y)) . V /(y) 



k - y\ l 



dy (l - 6> R (x - y)) e 



x — y 



■V, 



+ / dye" 



»k-«r / x — y 



\x - y| 2 



V/(y) 



x - y\ 2 v 
x-y 



x-y 



\x - y| 



x - y\ 2 
2 • V^(x-y)| 



V/(y) 



(6.20) 
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Therefore 



(IV) < CtU dx ( dy (1 - 6 R (x - y)) 



|V/(y)| 



\x - y 



+ Ct*[ dx[ dy (1 - 9 R (x - y)) 



|V 2 /(j/)l 
\x — y\ 2 



+ Cl-[ I dx [ J dy\V9 R (x-y)\ 
I -Or 



< C 9 , s t2 ||V/||_3s_ 
3+s 



+ Cfl.t* IIV/I 



^q,s <"■ 1 1 v j 1 1 3 S 

3+s 



|v/(y)l 

I 1 9 

\x — y\ A 

- Cn T ih II V 2 / 



(2 + A-IW 

^3^g s' 



(6.21) 



V g r ' 



1 

"g,s 771 3" 



|v/| 



i 

£5 



3+s 



<?,r" 3 3 



i?" 



|v 2 /| 



for 3/2^s^g^oo, g / 3/2 and 1 ^ r < 3^/(3 + q). Here Young's inequality and the scaling 
properties (|6.8p and (|6.9p have been used. 

Summarizing, 
||e lA '/|L^ (/) + (H) + (m) + (IV) 



< c, 



it « 



3 



+ 



1 + 3-3 



+ 



i 

£5 



+ C, 



1 

is 



3 _ 3 



|v 2 /| 



1 + 3_3 



i 

+ C <J,S TT3 3~ 



|v/| 



(6.22) 



for any i? > 0, 3/2 ^s^g^oo, q ^ 3/2, 1 ^ r < 3g/(3 + g). Optimizing the choice of R leads to 
R = Vt, so that flEZZD reads 



£aU g- 1 



+ ||V/||^)+- 3 

£3V T g^ 



T l|v 2 /|| r . 



(6.23) 



For the r.h.s. of (|6.23p to stay bounded in time, we need 1 ^ r ^ 2g+3> which requires q ^ 3. This 
completes the proof of the proposition. □ 



A Pointwise bounds on the two-body wave function 

In this section we investigate the boundedness properties in time and in space of the evolution of 
the two-body wave function (p® 2 when the two particles are coupled by the interaction Vn(x\ — X2), 
that is, where the dynamics is generated by the Hamiltonian 

f)£' 2) = -A x - A 2 + V N (x x - x 2 ) 

acting on L 2 (M 3 x R 3 ,dxidx2), as defined in (|2.ip . In spirit, these bounds are similar to the ones 
proved in Proposition 15.31 (in particular in (|5.28p for s = 00; the initial data however, is different). 
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Proposition A.l. Let V be a non-negative, smooth, spherically symmetric, and compactly supported 
potential. Let V N {x) = N 2 V(Nx) and f)^' 2) = -Ai - A 2 + V N { Xl - x 2 ). Let ip t = e-^'V 52 , for 
some (p G L 2 (R 3 ). Then, for every a > 3, there exists C > such that 

UtU < C||^||l )00)Q logiV (A.l) 

where 

4 

ibikoo,a = E IN QV "vllL- ( A - 2 ) 

m=0 



Proof. Let ip(r], x) = ip(r] + x/2, r] — x/2), and 

far], x) = ipt(v + x/2, J] - x/2) = (V iA " /2 *e-^*^) (77, x) 
where the operator f)jv = — 2A + Vn{x\ — x 2 ) only acts on the relative variable x. Then 

H^lli°°(R6,d:zid:E2) = ll^*lli°°(K 6 ,dr;da;) ^ C SU P H^*("> x ) Wm(M?,dn) 



2 r 

< CV sup / d V \V™Mv,x)\ 

2 _ (A.3) 



m=0 
2 



m=0 17 R 

It is useful to switch to the macroscopic coordinate X = Nx. With the short hand notation 

1>™(X) := 4 m) (§), with (A.4) 

^ m \x) := V™j( V ,x)=V™(v(ri + 1)^-1)), (A.5) 



Eq. (|A.3j) reads 



MlL°°(R 6 ,da;idiE2) 

m=0 



because f)jv = A r2 (— 2Ax + V(X)). To bound the r.h.s. of the last equation, we use the modified 
dispersive estimate (|6.4p with s = oo. We obtain 



II e-^^X^AX) < c (linvggllL + llVxnvSSlll + liv^nvjgll j) ' (A ' 7) 

for every m = 0, 1, 2, where f2 is the wave operator associated with —A + ^V(x) = f)/2, as defined 
in Proposition ID.ll 

To estimates the terms on the r.h.s. of ()A.7p we use the same strategy followed in the proof of 
Proposition 15.31 However, since the initial data considered here (the wave function VtJ™?!) * s different 
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from the one considered in Proposition 15.31 the bounds for the r.h.s. of (|A.7|) cannot be directly 
inferred from the analogous estimates for the terms of the r.h.s. of (|5.29p . 

By Yajima's bound (|D,2p . we have 



(m) I 



Analogously to the treatment of [|Vfl*u;^A||t> after (|5.36p . we obtain 

C 



llVx^Slls < _ /a (11 W$ 11^ + 11 W$[|a 



Since 



H^^SlI 3 = N*\\A x ^\\s 



Mi 



AT 113 



— HA ii>( m )|U 



(A.8) 

EG (0,1). (A.9) 

(A.10) 
(A.ll) 



one has 



T ( iV £ || A^ m) |U + ^ II A ^v m) h + 



,M|i 

T) II 00 



e(l-e)3 v ' ^ AT 

Similarly to the treatment of || V 2 0*u;^a||z in (|5.50j) . we find 

||VjQ><$|| f < c(||AxVSll| + II^Sll|) < c(||A^H|| | + ||^f) 

Using (EOl) . (|ATl2j> . and (|A"13"1) to bound the r.h.s. of (TO) , we find 

I -iljjV 2 t ./.(«*) ||3 



(A.12) 



(A.13) 



< C 



A £ 



,e(l-e) 



A,^ m) H 2 ^ + ||A^||| + A^H! + ||^ m) ||L) • 

(A.14) 



Now it remains to integrate the r.h.s. of the last equation over rj G R . By the assumptions on (p, it 
follows that 



|VV(s)| < IMk 



i 



(A.15) 



for any i/ = 0, . . . , 4 and some a > 3. This means that for any p € [1, +oo) and n,m G {0, 1, 2} one 
has 



l V "^ m) || p 



dx y?v- p U + - hy 



1 

p \ p 



m+n 
/i=0 



dx 



i 

p\p 



(A.16) 



< CNI<oo,<«) 



dx 



1 



(V + %) ap (v ~ f ) ap 
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and the above integral can be estimated as 



11 C 

dX (v+%) ap {r]~%) ap ^ (A ' 17) 



as long as ap > 3, which is always the case due to the assumption a > 3. Then 

II II 4 

\K4 m) \\l < C ilgpg n,m 6 {0,1,2} (A.18) 

independently of p G [1, +oo]. 

Due to (|XT8|) . when plugging ([O]) . (]XT2]) . and (fQ3]) into the r.h.s. of estimate (|X7j) . and the 
latter into the r.h.s. of (|A.3p . integrability in 77 is guaranteed by the assumption a > 3 and one gets 

Woo < HH,oo,a -r £€(0,1). (A.19) 

e(l-e)3 

After choosing e = (log A) -1 , estimate (jA.ip is proved. □ 

B Estimates for the energy of a factorized data 

In the following lemma we prove that the expectation of the Hamiltonian iTjy i n a factorized state 
(p® N , for some 99 G -ff 2 (M 3 ), is of the order A. This estimate is used in Section[3j Moreover, we also 
show that the expectation of (in the state ip® N ) is of the order A 3 . 

Lemma B.l. Let Hn be defined as in il.iy\) and let (p G -£f 2 (R 3 ). Then 

Inn \U® N ,H N V® N ) = \\^v\\l + WhMl V^G^(M 3 ), (B.l) 

^^(^,^0= ilMllrill V^Gtf 2 (M 3 ). (B.2) 



Proof. We will prove equation (|B.2|) . the proof of (|B.1|) is similar but simpler, so we will omit it. 
Using the permutation symmetry of ip® N , we find that 



A 3 



A(A - 1) 



< A -1 ((^® JV , AiA 2 ^ JV ) + A- 2 (^ 7V , A 2 <^ v ) 



iV a2,„®AT\ 



+ <<^ JV , (- Ax) V N (x 2 - x 3 )^ v ) + I <^ v , Ax VivCx! - x 2 )<p 
+ N(<p® N , V N ( Xl - x 2 )V N (x 3 - x A )<p® N ) 
+ (<p® N , V N ( Xl - x 2 )V N (x 2 - x 3 )<p® N ) . 

Since <p G H 2 (M. 3 ) and \\ip\\ = 1, it follows that 

(<p® N , AtAw®") < \\<p\\ 4 Hl < 00 and 

Moreover, using Lemma lB.21 we have 

{<p 9N , (-A 1 )V N (x 2 - x 3 )<p® N ) < C \\V N h \\ipf m < CN- 1 \\V\h \\tpf Hl 



J8>N\ 



1 (g>N A 2 ®iV\ ^ 11 nz , 

{<P , A^ ) < \\<p\\ H 2 < 00 



2 

H 2 



(B.3) 



(B.4) 
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and, with a Schwarz inequality and Sobolev embedding, 

(<p® N , A 1 V N (x 1 - x 2 )^ N ) | < (<p® N , A 1 V N (x 1 - x 2 )A w ® N } 1/2 &® N , V N ( Xl - x 2 )<ffi 

< CN~ l \\V\\i \\ipf H2 ■ 
Finally, using again Lemma lB.21 we observe that 

(tp® N ,v N ( Xl -x 2 )v N (x 3 -x 4 )<p® N ) <c\m 8 h1 \\Vn\\1 < cn- 2 \\v\\ 2 m\ 8 h1 

and that 

(<p® N , V N {x x - x 2 )V N (x 2 - x 3 )tp® N ) < C \\<pt Hl \\V\\ 2 3/2 . 
Inserting all these bounds in the r.h.s. of (|B.3p . it follows that 



iiV\l/2 



lim — 7 



(^ N ,H 2 N ^ N ) - ^-^{^ N ,V 2 ( Xl -x 2 )^ N ) 



(B.5) 



(B.6) 



Eq. (|B.2p now follows because 

-±^(^,Vi( Xl -x 2 )^ N ) = ± f dx 1 dx 2 N^V 2 {N{x 1 -x 2 ))\^x l )\ 2 \^x 2 )\ 2 ^^ 



I Il4 

ML 



as N —>■ oo (the convergence follows by a Poincare inequality, since <p £ H 2 (M. 3 )). □ 
Lemma B.2 (Sobolev-type inequalities). Let ip £ L 2 (R 6 , dxidx 2 ). If V £ L 3//2 (M 3 ), we have 

\^,V{ Xl -x 2 )^)\ < C||y|| 3 / 2 (1 - A x )^) . (B.7) 

IfV € L^IR 3 ), then 

K^Vfa-x^)] <C||F||i(^(l-Ai)(l-A 2 )V) (B.8) 

The first bound follows from a Holder inequality followed by a standard Sobolev inequality (in 
the variable x\, with fixed x 2 ). A proof of the second bound can be found, for example, in [1] [Lemma 
5.3]. 



C Properties of the one-body scattering solution 1 — uj{x) 

Lemma C.l. Let V : M 3 — > IR be non negative, smooth, spherically symmetric, compactly supported 
and with scattering length a. Let 1 — tu(x) be the solution of 



A+±V)(l-w) = 



with uj(x) — ► as \x\ — > oo . 



Then there exist constants C m ,C depending on the potential V such that 



\u(x)\ ^Conr 
I en | 

\v m oj(x)\ ^ a 



771 < x \m+l 



\uj(x)\ (C < 1) 

\v m cj(x)\ < c m 



(C.l) 

(C.2) 
(C.3) 



for every nonnegative integer m. As a consequence, for every p £ [l,oo] and every nonnegative 
integer m such that p(m + 1) > 3, one has 



|V m w||„ < C mt p < OO. 



(C.4) 
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Proof. Inequalities (|C.2p and (|C.3|) follow immediately from the fact that, out of the support of the 
potential V one has 

uj{x) = ^- \x\ > R, (C.5) 

while, inside the support, V m u; is bounded by elliptic regularity and compactness for any nonnegative 
m. The fact that the constant C in (|C.2|) is strictly smaller than 1 is proved in Lemma B.l of [2\. □ 

By scaling, one immediately has the following. 

Corollary C.2. Let V be as in Lemma \C.1[ Let Vn(x) = N 2 V(Nx) and let 1 — ujn be the corre- 
sponding solution of the zero-energy scattering equation. Then 

lo n (x) = lu(Nx) (C.6) 

whence 

\\V m io N \\ p = N m -l\\V m u;\\ p (C.7) 
under the same condition for the validity of fC-^| ), 

D Properties of the wave operator Q 

We denote by Q the wave operator associated with the one-particle Hamiltonian h := —A + \V , 
that is, h = |rj in our previous notation. Its existence and most important properties are stated in 
the following proposition. 

Proposition D.l. Suppose V > 0, with V G L X (IR 3 ). Then: 

i) (Existence of the wave operator). The limit 

Q = s- lim e iht e iAt 

t^oo 

exists. 

ii) (Completeness of the wave operator). Q is a unitary operator on L 2 (R 3 ) with 

n* = n- 1 = s - lim e-^e-M 

t— >oo 

Hi) (Intertwining relations). On D{h) = D{— A), we have 

n*hn = -A (D.i) 

iv) (Yajima's bounds). Suppose moreover that V(x) < C(x)~ a , for some a > 5. Then, for every 
1 <P < oo, O and Q* map L'P(R 3 ) into LP(R 3 ), that is, 



W^Wlp^lp < °° f or all 1 P ^ °° • (D-2) 
If moreover V 6 C k (M. 3 ), we have 

/or a// m < k. 

Proof. The proof of i), ii), and iii) can be found in [7]. Part iv) is proved in [101 E]. □ 
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E Correlation structure and Gross-Pitaevskii equation 



As remarked in the introduction, the correlation structure developed by the solution to the iV-particle 
Schrodinger equation VPjvt = e~ lHNt ^^ (with Hjq defined as in for initial states exhibiting 

complete Bose-Einstein condensation plays a very important role in the derivation of the Gross- 
Pitaevskii equation (|1.1U|) ; more precisely, the emergence of the scattering length in the coupling 
constant in front of the nonlinearity is a consequence of the presence of the short scale correlation 
structure. The goal of this appendix is to explain this connection in some more details (see also 
Section 3 in 121). 



From the Schrodinger equation (jl.3p . it is simple to obtain an evolution equation for the one- 
particle density 

idtiff t {xr,^i) = (-A Xl + A X A 7^ ) t (x i ;x / 1 ) 

f (2) (K1) 

+ (N-1) / dx 2 (V N (xi -x 2 ) - V N (x[ -x 2 )) T$ t (xi,x 2 ;x' 1 ,X2) ■ 

This is not a closed equation for 7^ because it also depends on the two-particle density 7^ as- 
sociated with *kN,t (the two-particle density is defined similarly to (jl.4p . integrating however only 
over the last (N — 2) particles); actually, (lE.lj) is the first equation of a hierarchy of N coupled 
equations, known as the BBGKY hierarchy, for the marginal densities of ^N,t- From 7&I — > \ift)((pt\ 
as N — > 00, it follows that jffi t — > \ip t )(ip t |® 2 . If we replace, in (|E.ip . the densities 7^ and j^' t 
by these limit points, and if we replace (N — I)Vn(x) — N 3 V(Nx) by its (formal) limit b5(x) with 
b = j V, we obtain a closed equation for the condensate wave function cpt, which has the same form 
as the Gross-Pitaevskii equation (jl.lOj) . but with a coupling constant in front of the nonlinearity 
given by b instead of 87ra. The reason why this naive argument leads to the wrong coupling constant 
is that the two-particle density 7^ contains a short scale correlation structure (inherited by the 
iV-particle wave function ^N,t) which varies on exactly the same length scale N^ 1 characterizing the 

interaction potential. Describing correlations by the solution to the zero energy scattering equation, 

(2) 

we can approximate, for large but finite N, the two-particle density ^ N t by 

^ t (x 1 ,X2;x' 1 ,x 2 ) ~ (1 - un(xi ~ - u} N (x[ - 4))^(^l)^(2 ; 2)^i(^lM(4) • 

Inserting this ansatz in the second term on the r.h.s. of (|E.l|h and using (jl.8p . we obtain the correct 
Gross-Pitaevskii equation for (ft. The emergence of the scattering length in the Gross-Pitaevskii 
equation is therefore a consequence of the singular correlation structure developed by *&N,t (and then 
inherited by the two particle marginal density). 
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